Cooperative Jahn- Teller Distortion in PrC>2 
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We report neutron diffraction data on single crystal Pr02 which reveal a cooperative Jahn- Teller 
distortion at Td = 120 ± 2 K. Below this temperature an internal distortion of the oxygen sublattice 
causes the unit cell of the crystallographic structure to become doubled along one crystal axis. We 
discuss several possible models for this structure. The antiferromagnetic structure below Tn = 13.5 K 
is found to consist of two components, one of which shares the same doubled unit cell as the distorted 
crystallographic structure. We also present measurements of the magnetic susceptibility, the specific 
heat capacity and the electrical conductivity of PrC>2. The susceptibility data show an anomaly at 
a temperature close to Td. From the specific heat capacity data we deduce that the ground state 
is doubly degenerate, consistent with a distortion of the cubic local symmetry. We discuss possible 
mechanisms for this. The conductivity shows an activated behaviour with an activation energy 
£ a = 0.262 ± 0.003 eV. 



PACS numbers: 61.66.Fn, 75.25,+z, 75.30.Cr, 75.40.Cx 
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I. INTRODUCTION 

In recent years there has been a resurgence of in- 
terest in phenomena associated with orbital degrees 
of freedom^ such as orbital ordering, orbital waves 
( "orbitons" )* and unusual spin-orbital liquid ground 
statesi&i Much of the work has been stimulated by stud- 
ies of the perovskite-based transition-metal oxides, e.g. 
the colossal magnetoresistance manganites, whose prop- 
erties can be strongly influenced by the extent to which 
orbital degrees of freedom are quenched by Jahn- Teller 
distortions of the lattice or coupled to other electronic 
degrees of freedom. 

Jahn- Teller and orbital phenomena also occur in com- 
pounds containing localized 4/ and 5/ electrons. Among 
the simplest of these, the nuorite-structure actinide diox- 
ides UO2 and Np02 have been investigated for many 
years and are now understood to exhibit complex or- 
dered phases at low temperatures involving coupled elec- 
tric and magnetic multipoles as well as (in the case of 
UO2) a lattice distortion. 

Unusual magnetic effects have also been found in the 
fluorite-structure lanthanide dioxide Pr02- Some years 
ago, Pr02 was found to have an anomalously small or- 
dered moment in the antiferromagnetic phased Very re- 
cently, we discovered a broad continuum in the magnetic 
excitation spectrum of Pr02 probed by neutron inelastic 
scattering, which we ascribed to Jahn- Teller fluctuations 
involving the orbitally-degenerate 4/ ground state and 
dynamic distortions of the latticed In a separate neu- 
tron diffraction experiment we found evidence that the 
antiferromagnetic structure contains a component with 
twice the periodicity of the accepted type-I magnetic 
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In this paper we follow up our previous work&S with 
further experimental results on the electronic and mag- 
netic behaviour of Pr02- Our main results come from 
neutron and x-ray diffraction measurements of the crys- 
tallographic and magnetic structure of Pr02 in the tem- 
perature range 2-300 K. These reveal the existence of 
an internal distortion of the fluorite structure below 
Td = 120 K and a related distortion of the antiferro- 
magnetic structure below Tn = 13.5 K. We also report 
measurements of the magnetic susceptibility, the specific 
heat capacity, and the electrical conductivity of Pr02, 
which allow us to determine the effective paramagnetic 
moment of the Pr ion, the degeneracy of the 4/ electron 
ground state, and the activation energy of the charge 
carriers. 



II. SAMPLE PREPARATION 

Measurements were made on single crystal and poly- 
crystalline samples of Pr02- The single crystals used for 
the neutron diffraction and conductivity measurements 
were taken from a batch prepared several years ago by a 
hydrothermal procedure. 1 — The sample used for neutron 
diffraction was the largest of the batch, with a mass of 
ss 1 mg, while the sample used for the conductivity mea- 
surements was somewhat smaller. Both were irregularly 
shaped. Polycrystalline samples were used for the x-ray 
diffraction, magnetic susceptibility and specific heat ca- 
pacity measurements. These were prepared by oxidation 
of commercially obtained PrgOn. The starting powder 
was first baked in air at 1000°C for ~ 11 hours, then 
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annealed in flowing oxygen at 280° C for 30 days. The 
powder was reground approximately once a week during 
the annealing. Synchrotron x-ray diffraction showed that 
the final product contained ~ 1% of residual PigOn. 



III. STRUCTURAL INVESTIGATION 

We investigated the crystallographic and magnetic 
structure of PrC>2 by performing a neutron diffraction 
study on a single crystal sample. The experiment was 
carried out on the DfO four-circle diffractometer at the 
fnstitut Laue-Langevin. We used an Eulerian cradle for 
crystal orientation and a position sensitive detector. The 
latter was kept in the scattering plane at all times. The 
diffractometer configuration was as follows: vertically 
curved pyrolytic graphite (002) or Cu (200) monochro- 
mator, no collimators, circular aperture of diameter 6 mm 
before the sample and rectangular aperture of dimensions 
20 mm x 25 mm before the detector. We used incident 
neutron wavelengths of 2.3575 A (Ai) and 1.2579 A (A2) 
from the pyrolytic graphite and Cu monochromators re- 
spectively, with a pyrolytic graphite filter for the former 
to reduce the half-wavelength contribution. The larger 
wavelength provided greater flux, while the smaller pro- 
vided greater coverage of reciprocal space. 

The single crystal sample was mounted on a thin alu- 
minium pin and aligned such that the [110] direction lay 
along the axis of the pin. The pin was attached to the 
Eulerian cradle inside a helium flow cryostat. 



A. Measurements 

The quality of the crystal was checked by performing 
crystal rotation scans (w-scans) of the strongest struc- 
tural Bragg reflections at room temperature. The peak 
widths were found to be within the experimental resolu- 
tion. 

To examine the crystallographic structure of Pr02, as 
many structural Bragg reflections as possible were mea- 
sured by w-scan at the accessible positions in reciprocal 
space satisfying the selection rules for the fluorite struc- 
ture (h, k and I all even or all odd). This set of reflections 
was measured at Ai at temperatures of 2K, 20 K, 90 K, 
150 K and 300 K. At some temperatures the full set was 
measured at both neutron wavelengths (Ai and A2). 

The intensities of the strongest reflections (e.g. (220), 
(400)) were found to change rapidly with temperature. 
The temperature dependence of these reflections was 
measured more carefully by performing w-scans at 2 K 
intervals between 2K and 130 K, then at 5K intervals 
between 130 K and 300 K. Figure ^ shows the integrated 
intensity of the scan as a function of temperature. A 
sharp rise in intensity can be seen below 120 K. We found 
that the rise was greatest for the strongest reflections and 
more pronounced at the larger wavelength Ai, suggesting 
a sudden reduction in extinction. A possible explanation 



for this observation is a reduction in grain size caused by 
the strain associated with a structural distortion. 
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FIG. 1: Temperature dependence of the (220) reflection. 



The temperature dependence of the lattice parame- 
ter, estimated from the variation of the Bragg angle with 
temperature for the (400) reflection, also shows a sharp 
change at ~ 120 K (see Fig.0 closed circles). Although 
the absolute value of the lattice parameter obtained by 
this method is inaccurate, as it relies on the crystal be- 
ing perfectly aligned, the data nevertheless give a good 
qualitative indication of the way the lattice parameter 
changes with temperature. It is seen to decrease linearly 
from 300 K to ~ 120 K, below which it becomes approx- 
imately constant, dipping slightly at ~ 18 K. The open 
circles in FigureElshow data taken from a high-resolution 
x-ray powder diffraction study on PrC>2 performed at 
the Swiss-Norwegian Beamline BM1B at the European 
Synchrotron Radiation Facility (ESRF). There are fewer 
points in this plot, but the absolute value of the lattice 
parameter is much more accurate because it has been ob- 
tained by Rietveld refinement of the diffraction pattern. 
The temperature variation found by x-ray diffraction is 
consistent with that found by neutron diffraction. 

To probe the magnetic structure, as many magnetic 
Bragg reflections as possible were measured by w-scan 
at the accessible positions in reciprocal space satisfying 
the selection rules for the antiferromagnetic (AFM) type- 
I magnetic structure (h, k and I a mixture of even and 
odd). This set of reflections was measured at tempera- 
tures of 2K and 20 K (below and above Tjj) using only 
the larger wavelength Ai , since this provided higher flux. 

In addition to the above measurements we also per- 
formed cj-scans at a number of positions in reciprocal 
space with half-integer Miller indices. The object of these 
measurements was to search for peaks seen previously in 
a powder diffraction experiment £ The scans were per- 
formed at temperatures of 2K, 20 K and 150 K. 

At T = 2 K, we found reflections at positions satisfying 
the selection rule h = n + 1/2, k = odd, I = even, where 
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FIG. 2: Temperature dependence of the lattice parameter in 
Pr0 2 . 



n, k and I are integers, and h, k and I can be commuted. 
At T = 20 K (above Tn), we found reflections at exactly 
the same positions, except that those with I = were 
no longer present. The strongest half-integer reflection 
we observed was ~1% of the intensity of the strongest 
fluorite reflection. At T — 150 K, all the half-integer 
reflections had disappeared. Figure [3] shows w-scans of 
the (| 10) and (±14) reflect ions. These are typical of the 
two observed sets of half-integer reflections: those with 
I = and those with I ^ 0. 
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FIG. 3: tj-scans of typical half-integer reflections (detector 
counts are normalised to a fixed incident beam monitor count 
of 1 x 10 6 , corresponding to a counting time of ~ 100 seconds). 
Those reflections with 1 = are absent above Tn. Those with 
I 7^ are still present above Tn. All are absent at T = 150 K. 



We measured the temperature dependence of the (^10) 
reflection by counting at the peak centre at a series of 
temperatures between 2 K and 20 K. The temperature de- 
pendence of the stronger (^14) reflection was measured 
by performing w-scans at a series of temperatures be- 
tween 2K and 150 K. 

Figure ^compares the temperature dependence of the 
(±10) reflection with that of the magnetic reflection 
(100). To within experimental uncertainty, both reflec- 
tions disappear at the antiferromagnetic ordering tem- 



perature Tn- We were able to determine Tn for both re- 
flections by fitting the data using order parameter curves 
of the form 



I oc (T N - T) 



(1) 



where I is the peak count and j3 is the critical exponent. 
For the (100) reflection we found T N = 13.41 ± 0.04 K, 
while for the (±10) reflection we found T N = 13.5±0.2K. 
Since these are in excellent agreement, it is likely that 
both reflections are magnetic in origin and belong to the 
same magnetic phase. The existence of half-integer mag- 
netic Bragg reflections indicates that the magnetic unit 
cell is doubled along one crystal axis. 
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FIG. 4: Temperature dependence of the (|l0) and (100) in- 
tensities (detector counts are normalised to a fixed incident 
beam monitor count of 4 x 10 6 , corresponding to a count- 
ing time of ~ 8 minutes in the upper graph and ~ 6.5 min- 
utes in the lower graph) . Order parameter curves of the form 
I oc (Tn — T) 2/3 have been fitted to the data to determine 
Tn. Note that the temperature dependence of the (100) re- 
flection was actually measured during an earlier experiment 
performed on the same instrument with the same crystal, but 
with a different sample environment, so the intensities in the 
two plots should not be compared directly. 



Figure shows the temperature-dependence of the 
(±T4) reflection. It is found to disappear at Td = 
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120 ± 2K, which is the same as the temperature below 
which we observed a sharp rise in the intensities of the 
large structural reflections. This strongly suggests that 
the (±14) reflect ion is structural in origin. The exis- 
tence of half-integer structural reflections indicates that 
the crystallographic unit cell is also doubled along one 
crystal axis. 
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FIG. 5: Temperature dependence of the ( i 14) intensity. An 
order parameter curve of the form I oc (Td — T) 2/3 has been 
fitted to the data to determine Td. 



A re-examination of previous data taken with a pow- 
der sample on the POLARIS diffractometer at the ISIS 
facility^ has now also revealed half-integer structural re- 
flections. Although two of the magnetic half-integer re- 
flections had been seen previously, the structural reflec- 
tions had not been noticed, due to the nearby presence of 
other peaks of similar intensity, which were probably due 
to scattering from the sample environment. The x-ray 
powder diffraction study mentioned above also confirms 
the presence of half-integer structural reflections. 



B. Structure analysis 

To identify the nature of the crystallographic struc- 
ture below Td = 120 K and the magnetic structure be- 
low Tn = 13.5 K, we performed a detailed analysis of the 
intensities of the observed integer and half-integer struc- 
tural and magnetic reflections at low temperatures. 

After background subtraction, each aj-scan was in- 
tegrated using the extended trapezoidal ruleii (a sim- 
ple method of numerical integration), corrected for the 
Lorentz factor, normalised to the monitor count and 
the intensity of the (311) reflection, and multiplied by a 
wavelength-dependent scale factor to enable direct com- 
parison between the intensities of the observed struc- 
tural reflections and their calculated structure factors. 
The (311) reflection was used for normalisation because 
it was one of the weakest reflections accessible at both 
wavelengths, so it was little affected by extinction. In 



addition, its structure factor contained no contribution 
from the oxygen ions, so it was unaffected by the distor- 
tion. Finally, the corrected and scaled intensities were 
averaged over symmetry-equivalent reflections to guard 
against the possibility that a chance uneven distribution 
in symmetry-equivalent structural or magnetic domain 
populations might affect the measured intensity of any 
given reflection (although this was unlikely since the ma- 
jority of the symmetry-equivalent reflections were equal 
in intensity to within experimental uncertainty). 

We then considered distortions of the fluorite crystallo- 
graphic structure and the type-I antiferromagnetic struc- 
ture which would give rise to structure factors in agree- 
ment with the relative intensities and selection rules of 
the observed reflections. Such structures were identified 
by trial and error. However, in the case of the crystal- 
lographic structure we were guided by the results of the 
x-ray diffraction study, which revealed no splitting of the 
reflections caused by the fluorite structure. This indi- 
cated that the Pr lattice remained undisturbed, and any 
distortion must therefore be due to an internal rearrange- 
ment of the oxygen atoms. 

Once a suitable distorted crystallographic structure 
was found, we refined the overall magnitude of the oxy- 
gen displacements by a least squares method to achieve 
the best possible agreement between the structure factors 
and peak intensities. For the magnetic structure we used 
the ratio of the peak intensity to the magnetic structure 
factor for each reflection to obtain a value for the ordered 
moment /i of the Pr ion. 

The structure factors of the structural reflections are 
given by 



I^n(Q)| 2 = 



E 
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(2) 



where the summation index j runs over all the atoms in 
the unit cell, bj is the nuclear scattering length of the j'th 
atom averaged over all of its isotopes, Vj is the position 
of the jth atom within the unit cell, and e~ w j\Q' T > is 
the Debye- Waller factor (we set this equal to 1, since the 
measurements were made at low temperatures) and () 
denotes an average over all symmetry-equivalent struc- 
tural domains. The scaled, integrated intensities of the 
structural reflections are equal to |Fn(Q)| 2 . The struc- 
ture factors of the magnetic reflections are given by 



I^m(Q)| 2 = ((<W - Q a Qp) i^(Q)<*(Q) 

a/3 



(3) 



where the summation indices a and (3 run over the carte- 
sian co-ordinates x,y and z, 5 a p is the Kronecker delta, 
* denotes the complex conjugate, Q a is the a-component 
of the unit scattering vector and () denotes an average 
over all symmetry-equivalent magnetic domains. In the 
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dipole approximation F^(Q) is given by 



*m(Q) = /(Q)£ 



<e iQl V 



-Wj(Q,T) 



(4) 



where the summation index j runs over all the magnetic 
atoms in the magnetic unit cell, /t" is the ct-component 
of a unit vector in the direction of the magnetic moment 
of the jth magnetic atom, Vj is the position of the jth 
magnetic atom within the magnetic unit cell, /(Q) is the 
magnetic form factor of the Pr 4+ ion and e~ Wj ^ ,T ^ is the 
Debye- Waller factor (set equal to 1). The scaled, inte- 
grated intensities of the magnetic reflections are propor- 
tional to (2lo) 2 (J^) 2 |f M (Q)|2 j where (Tio) 2 = 72 .4mb 

(= 7.24 fm 2 ) and /j,b is the Bohr magneton. 

For reference we note the structure factors of the fluo- 
rite structure here: 



336 fm , h, k, I all odd 
|F N (Q)| 2 = { 791 fm 2 , h + k + l = 2n 
4194fm 2 , h + k + l = 4n. 



(5) 



These were computed using the fluorite unit cell, which 
contains four Pr ions and eight O ions. For the half- 
integer reflections (both structural and magnetic) we 
used a unit cell that was doubled along one axis, thus 
containing eight Pr ions and sixteen O ions. In order to 
allow direct comparison with the structure factors of the 
fluorite structure we divided the structure factors Fn(Q) 
and Fm(Q) of the half-integer reflections, calculated in 
the doubled cell, by 2. 



Distorted crystallographic structure 



A full crystal structure refinement using the Cambridge 
Crystallography Subroutine Library (CCSL) yielded d = 
0.0664(17) A. 
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FIG. 6: Sheared structure, (a) The large pale spheres are 
praseodymium ions and the small dark spheres are oxygen 
ions. The oxygen cubes in the two halves of the doubled unit 
cell are sheared in opposite senses, but both have a shear- 
ing vector that is perpendicular to the direction along which 
the cell is doubled. The displacements of the oxygen ions 
have been exaggerated for clarity, (b) Scale diagram showing 
the distorted structure in plan view. The white circles are 
praseodymium ions and the black circles are oxygen ions. 



One of the simplest structures that possesses struc- 
ture factors in agreement with the intensities of the half- 
integer structural reflections is shown in Fig. |5{a). The 
Pr lattice is undisturbed, but the oxygen ions are each 
displaced by an amount d in a direction perpendicular 
to that along which the unit cell is doubled (in Fig. 
the doubling is along the x-direction and the displace- 
ments are along the y-direction) . The oxygen cubes in 
the two halves of the doubled unit cell are sheared in 
opposite senses, so we will refer to this structure as the 
"sheared" structure. The space group is Imcb (a variant 
of Ibam, space group 72), with the origin shifted by \ of 
the unit cell length along the orthorhombic a-direction 
(the x-direction in Fig.EJ. 

A simple least squares refinement of the oxygen 
displacement d, based on comparison of the domain- 
averaged structure factorsSi and the observed intensities 
of the half-integer reflections, yielded d = 0.0726(54) A. 
The structure factors and observed intensities at T = 
20 K are listed in Table [I] and a plan view of the struc- 
ture with displacements to scale is shown in Fig. H^b). 



An identical shearing of the oxygen cube was originally 
proposed for an internal distortion observed in UO2, 12 
although in that compound there was no evidence for a 
doubling of the unit cell. In UO2 the oxygen ions are dis- 
placed by 0.014 A from their fluorite structure sites, so 
the oxygen displacements in Pr02 are five times as large, 
and this is reflected in the high value of the transition 
temperature. It is now accepted that the oxygen con- 
figuration in UO2 forms a triple-q structured However, 
there is currently no reason to assume that the same is 
true in Pr02- 

It is interesting to note that there are two distinct Pr 
sites in the sheared structure, each of which occurs with 
equal frequency. The two sites have different surround- 
ing oxygen configurations, which are shown in Fig. [3 At 
one site the oxygens form a parallelohcdron (a polyhe- 
dron with two square faces, two rectangluar faces and 
two parallelogram faces), while at the other they form a 
polyhedron with four parallelogram faces and four trian- 
gular faces. 
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Reflection 


|IF N (Q)| 2 (fm 2 ) 


Intensity (fm 2 ) 


(3 10) 


0.0 


0.1 ± 0.1 


(fio) 


0.0 


0.3 ± 0.1 


(fio) 


0.0 


0.5 ± 0.2 


(312) 


5.1 


5.5 ± 0.3 


(§12) 


5.1 


5.0 ± 0.3 


(|12) 


5.1 


4.8 ± 0.5 


(112) 


5.1 


4.8 ± 0.6 


(2 14) 


19.9 


16.6 ± 0.9 


(§14) 


19.9 


20.5 ± 1.2 


(3 16) 


42.7 


42.3 ± 2.3 


(1 30) 


0.0 


0.4 ± 0.2 


(§32) 


5.1 


5.2 ± 0.6 


(§32) 


5.1 


7.1 ± 0.6 


(§32) 


5.1 


7.1 ± 1.0 


(132) 


5.1 


6.3 ± 1.1 


(134) 


19.9 


17.7 ± 1.5 


(§34) 


19.9 


22.0 ± 1.5 


(§36) 


42.7 


45.6 ± 3.0 


(152) 


5.1 


6.9 ± 1.2 


(§52) 


5.1 


8.7 ± 1.6 


(152) 


5.1 


4.9 ± 1.8 


(1^2) 


5.1 


12.4 ± 2.2 



TABLE I: Comparison between ||Fn(Q)| for the sheared 
structure and the intensities of the observed half-integer struc- 
tural reflections at T = 20 K. The Fn(Q) have been divided 
by 2, as described in Sec. 1111 Bl to allow direct comparison 
between the two columns. 
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FIG. 7: The two different Pr sites that occur in the 
sheared structure shown in Fig. |S| The white circles are 
praseodymium ions and the black circles are oxygen ions, 
(i) Parallelohedron. (ii) Polyhedron with four parallelogram 
faces and four triangular faces. 



consistent with the intensities of the observed reflections, 
it is not the only such structure. For example, a super- 
position of two structures identical to the sheared struc- 
ture, but with oxygen displacements of 0.0511 A in mu- 
tually perpendicular directions, gives similar agreement 
between the calculated structure factors and observed in- 
tensities. The overall displacement of each oxygen ion in 
that case is V% x 0.0511 = 0.0722 A. 

The "chirai" structure shown in Figgis also consistent 
with the intensities of the observed reflections, and has 
the additional advantages of having a lower Jahn- Teller 
energyM and a single Pr site. A chirai model is adopted 
for the oxygen displacement vector, but the structure fac- 
tors and the magnitude of the oxygen displacement are 
almost identical to that of the sheared structure (a sim- 
ple least squares refinement yielded d = 0.0726(54) A, 
whereas a full crystal structure refinement using CCSL 
yielded d = 0.0663(16) A). The space group is I4(l)/acd 
(space group 142), with Pr at positions 8a (0, |, |) and 
O at positions 16e (j + d, 0, |), where the origin is cho- 
sen to be at the centre of inversion, and c is the unique 
axis. 
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FIG. 8: Chirai distorted structure, (a) The large pale spheres 
are praseodymium ions and the small dark spheres are oxy- 
gen ions. The displacements of the oxygen ions have been 
exaggerated for clarity, (b) Scale diagram showing the chirai 
structure in plan view. The white circles are praseodymium 
ions and the black circles are oxygen ions. The + and — sym- 
bols indicate positive and negative displacements along the 
z-direction. 



Although the distorted structure described above is 
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2. Doubled magnetic structure 

The presence of half-integer magnetic reflections along- 
side stronger integer magnetic reflections indicates that 
the magnetic structure consists of two components: a 
primary component with the same unit cell as the flu- 
orite structure and a secondary component with a unit 
cell that is doubled along one crystal axis. This picture 
fits well with the model of the distorted crystallographic 
structure, in which the Pr lattice remains unchanged, 
but the oxygen sublattice undergoes an internal distor- 
tion which gives rise to a component of the structure with 
a doubled unit cell. 

The easiest way to analyse the magnetic structure is 
to consider the two components separately, so that the 
overall structure can be visualised by performing a vector 
addition of the two components of the magnetic moment 
for each Pr ion in the doubled unit cell. 

The primary component of the magnetic structure is 
the AFM type-I structure, shown in Fig. OH This gives 
rise to the integer magnetic reflections. Under ambient 
conditions it is not possible to distinguish between the 
three multi-q structures shown in the figure, but in view 
of the fact that the structural distortion doubles the unit 
cell along one crystal axis only, the single-q structure 
seems the most likely candidate. 
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FIG. 9: Type-I primary component of the magnetic structure. 
Neutron diffraction measurements under ambient conditions 
cannot distinguish between the multi-q structures shown, be- 
cause each gives rise to identical magnetic structure factors 
when averaged over symmetry-equivalent magnetic domains. 
However, since the crystallographic unit cell is doubled along 
one crystal axis, the single-q structure seems the most likely 
candidate. 



The configuration of the magnetic moments in the sec- 
ondary component of the magnetic structure was de- 
duced using a technique similar to that for the structural 
distortion. In our search for possible configurations we 
found it impossible to conceive of a magnetic structure 
that would give rise to reflections at I = alone. Sym- 
metry requires that if a structure gives rise to reflections 
at I = it will also give rise to reflections at all positions 
with I = even. This means that many of the magnetic 
half-integer reflections coincide with the structural ones. 
However, since little change was observed in the intensi- 
ties of the structural half-integer reflections below Tn (see 
Fig. EJl, it can be assumed that the magnetic intensities at 



these positions are much smaller than the structural in- 
tensities. We have therefore made no attempt to analyse 
the magnetic contribution at these positions. Our analy- 
sis of the secondary component of the magnetic structure 
is based solely on the half-integer reflections observed at 
1 = 0. 

The number of possible magnetic structures consistent 
with the observed half-integer intensities is quite large. 
Two of the simplest possible structures are shown in Fig. 
ITU1 Both possess identical magnetic structure factors, 
shown in Table [n] Figure HOf a) corresponds to the chi- 
ral model for the structural distortion shown in Fig. [SJ 
whereas Fig. llOf b) corresponds to the model shown in 
Fig. [S] The moments of the Pr ions in both magnetic 
structures point along directions perpendicular to the 
direction along which the unit cell is doubled. Similar 
structures in which some or all of the moments point 
along the doubling direction give poorer agreement with 
the relative intensities of the observed reflections, ft 
should be mentioned that the (§10) reflection had a high 
and rather sloping background, due to its proximity in 
reciprocal space to an aluminium powder line, so its in- 
tensity is somewhat unreliable. 
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FIG. 10: Possibilities for the secondary component of the 
magnetic structure, (a) All moments point along the y-axis. 
This structure is associated with the chiral model for the 
structural distortion shown in Fig. |H| (b) Half the moments 
point along the i/-axis and half along the z-axis. This struc- 
ture is associated with the model shown in Fig. |S| The mag- 
netic structure factors of both structures are identical after 
domain averaging. 



The magnetic moment of the Pr ion can be calculated 
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Reflection 


|^m(Q)| 2 


Intensity (fm 2 ) 


(s 10) 


1.5 


1.4 ± 0.1 


(fio) 


1.9 


1.7 ± 0.1 


(|io) 


1.8 


0.9 ± 0.3 


(|30) 


0.9 


1.1 ± 0.2 


(§30) 


1.0 


1.1 ± 0.2 



TABLE II: For the doubled component of the magnetic struc- 
ture the magnetic structure factors | i_Fivr(Q)| 2 (dimension- 
less) are compared with the magnetic intensities measured at 
T = 2 K. The Fm(Q) have been divided by 2, as described in 
Sec. 1111 Bl to maintain consistency and with the structural 
calculations. The ratio of magnetic intensity to magnetic 
structure factor is constant, and can be used to determine 
the magnetic moment of the Pr ion as described below. 



for each of the two components of the magnetic structure 
as follows: 



JL 

Mb 



/m(Q)| j f i n(Q') 



(3?) /n(Q')I^m(Q)| 2 



(6) 



where fj, is the magnetic moment of the Pr ion in one 
component of the magnetic structure, Im(Q) is the inte- 
grated intensity of a magnetic reflection with reciprocal 
lattice vector Q, |Fm(Q)| 2 is the magnetic structure fac- 
tor of this reflection, /n(Q') is the integrated intensity of 
a nuclear reflection with reciprocal lattice vector Q' and 
|-Pn(Q')| 2 i s the nuclear structure factor of this reflection. 

For each of the magnetic reflections we calculated the 
magnetic moment of the Pr ion from the observed mag- 
netic intensity and the intensity of the (311) nuclear re- 
flection. The resulting values of the Pr moment were av- 
eraged to reduce the uncertainty. For the primary (AFM 
type-I) component of the magnetic structure we obtained 
/ii = (0.65 ± 0.02)^b- This is in good agreement with a 
previous value of (0.68 ± 0.07)/xb^ obtained by neutron 
diffraction from the same single crystal sample. Previ- 
ous powder diffraction experiments have yielded values 
of (0.6 ± O.I^b 7 and (0.572 ± 0.012)/x B ^ but it is likely 
that these are in error, since in both cases the Rietveld 
refinement assumed the crystallographic structure to be 
fluorite at low temperatures, and therefore took no ac- 
count of the structural distortion. For the secondary 
(doubled) component of the magnetic structure we ob- 
tained [12 = (0.35 ± 0.04)^b (we disregarded the unre- 
liable (| 10) reflection when averaging the values of ^2 
obtained from the half-integer reflections). The uncer- 
tainties on our values of fi\ and [12 are dominated by the 
uncertainty on the magnetic form factor of Pr 4+ . 



IV. BULK PROPERTIES 

We now present measurements of some of the bulk 
properties of PrC^: magnetic susceptibility, specific heat 
capacity and electrical conductivity. These provide sup- 
port for the findings of the neutron diffraction experi- 
ments described above, and reveal clues as to the origin 
of the structural distortion. 



A. Magnetic susceptibility 

We measured the magnetic susceptibility of PrC>2 us- 
ing a commercial SQUID magnetometer, with a powder 
sample of mass 270 mg. The measurements were made 
using the reciprocating sample option (RSO), with an 
applied field of H = 1 T. Data were taken while cooling 
in steps from T — 350 K to T = 2K, with a delay to 
allow temperature equilibriation at each step. A plot of 
the molar susceptibility (per mole Pr) is shown in Fig. II II 
The inset shows the inverse molar susceptibility. Three 
features are evident. First, there is a peak at T = 14 K, 
due to antiferromagnetic ordering. Second, there is a 
change in gradient at T — 122 ± 2K which is seen most 
clearly in the inverse susceptibility. Following our struc- 
tural investigation we can identify this with the internal 
distortion at Tb- This anomaly had not been noticed 
in previous measurements^*^ due to the intrinsic noise 
present. Third, there is an upturn in the susceptibility 
below T w 6K which is due to the presence of a small 
amount of PrgOn in the sample. We demonstrate this 
in Figure El which shows a plot of the magnetic suscep- 
tibility of PrgOn measured by AC susceptometry with a 
powder sample of mass 327 mg. Between 1.5 K and 2.5 K 
there is evidence for a magnetic transition with glassy 
characteristics. At 2 K the molar susceptibility (per mole 
Pr) of P^On is ~ 100 times that of PrC>2, so the up- 
turn below 6 K in the PrC>2 susceptibility is consistent 
with the 1% PrgOn impurity known to be present in the 
sample. There is no feature in the PrgOn susceptibility 
at Td- 

We used our measurement of the Pr02 susceptibility to 
calculate the effective paramagnetic moment of the Pr ion 
at various temperatures, assuming that the susceptibility 
obeys the Curie- Weiss law for an antiferromagnet above 



X 



30/c B (r + 6l)' 



(7) 



where ^ mo1 is the molar magnetic susceptibility, 9 is the 
Weiss constant, Na is Avogadro's number, /cb is Boltz- 
mann's constant and n e ff is the effective paramagnetic 
moment of the Pr ion. x mo1 is in c g s units (emu [mol 
Pr] -1 ), as dictated by convention, whereas the quanti- 
ties on the right hand side of Equation (JJJ are in SI units. 
The conversion factor from SI to cgs units is incorporated 
into the right hand side. 
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FIG. 11: Magnetic susceptibility of PrCh powder. Molar sus- 
ceptibility is plotted in the main graph and inverse molar 
susceptibility is shown in the inset. 
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FIG. 12: AC magnetic susceptibility of PreOn powder. The 
real and imaginary parts of the molar susceptibility are plot- 
ted in the main graph, and the inverse of the real part of the 
molar susceptibility is shown in the inset. The signal in the 
imaginary part below T = 2.5 K indicates dissipation accom- 
panying the transition. 



Over the temperature range T = 250-350 K we ob- 
tained fi e ff = 2.32/xb. This is close to the value of 2.54^b 
expected for a free Pr 4+ ion. Below Td the gradient of the 
inverse susceptibility first decreases then increases, caus- 
ing /j, e ff to decrease. This indicates a reduction in the 
number of magnetic degrees of freedom, which is prob- 
ably caused by a lifting of the degeneracy of the crystal 
field ground state by the structural distortion (it is gen- 
erally accepted that the ground state is a quartet in 
the fluorite phase) 

Although an anomaly at T = 122 K had not been no- 
ticed in previous measurements, the curvature of the sus- 
ceptibility trace below this temperature had been noticed 
previously by Kern;: 16 When attempting a crystal field 



analysis he found that his data could not be fitted well 
with a crystal field of cubic symmetry. He proposed that 
the oxygen ions surrounding the Pr ion did not in fact 
form a cube, but that the two at opposite ends of the 
body diagonal were displaced outwards along the (111) 
direction, causing a splitting of the T$ crystal field level 
into two doublets. Using this modified structure he ob- 
tained a good fit to his data over the whole temperature 
range 0-300 K. He calculated the splitting of the T$ level 
to be 28.8 meV. 

The significance of this result was later dismissed 15 
because the sample used did not exhibit a clear anti- 
ferromagnetic transition. Although a discontinuity was 
observed at 14 K, the susceptibility continued to rise at 
lower temperatures. The sample had been prepared from 
a starting material of PrgOn by annealing at a tem- 
perature of 360° C under 5 atmospheres of oxygen pres- 
sure. These conditions lie close to the border, on a 
temperature-pressure diagram, of the regions in which 
the PrC>2 and PreOn phases are stable, so it is likely 
that the sample was not single phase PrC>2. However, 
x-ray diffraction showed that it had a lattice parameter 
of 5.393 ± 0.007 A, which agrees well with the accepted 
value for PrC>2 of 5.393 ± 0.001 A, and is much less than 
the value for Pr 6 On of 5.468 A H This suggests that the 
sample was mostly PrC>2, but contained enough of the 
PrgOn phase to obscure partially the antifcrromagnetic 
transition. By adding proportions of our susceptibility 
data for PrgOn and PrC>2, and comparing the results 
with the data obtained by Kern, we estimate that his 
sample contained ~ 20% Pr 6 On. 

Since we believe that Kern's sample was predominantly 
PrC>2, and we have now obtained direct evidence for an 
internal distortion of the oxygen sublattice below Td = 
120 K, we believe that his original crystal field analysis, 
in particular the 28.8 meV splitting of the Tg crystal field 
ground state is still of some relevance below Td- 



B. Specific heat capacity 

We measured the specific heat capacity of PrC>2 by the 
relaxation method^ using a laboratory-built calorime- 
ter. The sample was a disc-shaped pressed pellet of 
PrC>2 powder of diameter 10 mm, thickness 1mm and 
mass 83 mg, which was mounted on the sapphire sam- 
ple platform of the calorimeter using a small amount of 
Apiezon grease (high thermal conductivity grease). 

The measurements were made over a temperature 
range 2.4 K to 23 K. We measured the heat capacity of 
the sapphire platform first, then the heat capacity of the 
sample and platform combined. We obtained the heat 
capacity of the sample alone by subtracting the former 
from the latter. Finally, we obtained the specific heat 
capacity of the sample by dividing by the sample mass. 

Figure IT31 shows a plot of the specific heat capacity of 
PrC>2 versus temperature. The experimental uncertainty 
is estimated to be ~ 5%. A lambda anomaly due to the 
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antiferromagnetic ordering is observed at ~ 13.5 K. This 
is superimposed on the contribution to the specific heat 
from vibrations of the crystal lattice. 

The specific heat capacity of a pressed pellet of CeC>2 
powder (commercially obtained) was also measured. 
This sample had mass 81 mg and was of similar shape 
to the PrC>2 sample. Since Ce02 is non-magnetic, but 
has the same crystal structure, a similar lattice parame- 
ter and similar formula mass to PrC>2, it provides a good 
estimate of the contribution of the crystal lattice to the 
PrC>2 specific heat. The CeC>2 data is plotted with the 
Pr0 2 data in Fig. El 



estimated using the CeC>2 data: 

f 23 c(CeQ 2 ) 17 ™-i v -i , nf rt 

Asiattico = / = dT = 1.7 JK kg . (10) 

Jo 1 

By subtracting (|1U|I from J5J we obtained the change in 
magnetic entropy: 

As magnct ic = As to tai - As latticc = 31.2 JK _1 kg _1 . (11) 
From Boltzmann's law, 



PrO, 



CeO, 



_ln 2_ 
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Temperature (K) 

FIG. 13: Specific heat capacity of PrC>2 and CeCh- The white 
circles are the Pr02 data and the black circles are the Ce02 
data. The Ce02 data provide a good estimate of the phonon 
contribution to the PrC>2 specific heat. 



= nfce ln g, 



(12) 



where n is the number of magnetic ions per unit mass 
(for PrC>2, n = 3.48 x 10 24 kg _1 ) and g is the degen- 
eracy of the ground state above Tn (below Tn the de- 
generacy is lifted). The entropy changes predicted for 
a quartet (<? = 4) and a doublet (g — 2) ground state 
are 66.7JK _1 kg _1 and 33.3 JK _1 kg _1 respectively. The 
measured change (|11|) suggests that the ground state is 
a doublet. This finding supports a distortion of the crys- 
tal structure, which would lift the degeneracy of the T 8 
quartet below Tb- An extension of the specific heat ca- 
pacity measurement up to room temperature, to look for 
an anomaly near 120 K, would be eminently worthwhile. 



C. Electrical conductivity 



To determine the degeneracy of the crystal field ground 
state in PrC>2, we used our measurement of the specific 
heat capacity to calculate the change in magnetic entropy 
on passing through the magnetic transition. The specific 
entropy is related to the specific heat capacity by the 
formula 

s(T)=J T ^dT, (8) 

where c is the specific heat and s is the change in spe- 
cific entropy between T = and a finite temperature T. 
The total specific entropy of PrC>2 (magnetic entropy + 
lattice contribution) was obtained from the specific heat 
data by dividing by T, then integrating by the extended 
trapezoidal rule>ii This gave 

A Stotal = r dT = 32.9 JK^kg- 1 . (9) 

Jo 1 

To obtain the change in magnetic entropy it was neces- 
sary to subtract off the lattice contribution. This was 



We measured the electrical conductivity of PrC>2 as a 
function of temperature using a simple probe designed 
for use with a helium cryostat. Two thin gold wires were 
attached to opposite sides of the single crystal sample 
using silver conducting paint (attempts to attach four 
wires to the tiny crystal were unsuccessful, due to the 
difficulty in achieving good electrical contact while keep- 
ing the contacts separate). The crystal was glued to the 
copper base of the probe with GE varnish to ensure good 
thermal contact. 

Four-terminal measurements were made of the resis- 
tance of the sample + gold wires (this was justified, since 
the resistance of the crystal was much greater than the re- 
sistance of the gold wires) over a temperature range from 
5.8 K to 288 K. The temperature was increased in steps, 
allowing the sample to equilibriate at each new setpoint. 
Figure E] shows a plot of current against temperature at 
a voltage of 3 V. Below T ~ 180 K the current became too 
small to measure reliably. We were unable to determine 
the absolute value of the electrical conductivity because 
the shape of the crystal was irregular and its dimensions 
were not known exactly. 

The inset to Fig.ll4lshows that the conductivity has an 
activated behaviour, at least over the temperature range 
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FIG. 14: Current I through a single crystal of PrC>2 at a 
constant voltage of 3 V. The conductivity is proportional to 
I. The inset shows In / against 1/T. 



180-288 K. To determine the activation energy we fitted 
our data using the Arrhenius form: 

' = 'oexp(-^), (13) 

where / is the current and E a is the activation en- 
ergy. We extracted E a from several runs, each mea- 
sured at a different constant voltage (between — 5V 
and +5V), and took the mean of these values to ob- 
tain E a = 0.262 ± 0.003 eV for the temperature range 
180-288 K. This is much smaller than the band gap of 
~ 2.5 eV predicted by band structure calculations^ and 
is also smaller than the activation energies of other PrOa; 
materials, which have been reported in the range 0.4- 

Possible mechanisms for electrical conductivity in 
PrC>2 are (i) ionic conductivity, (ii) intrinsic electronic 
conductivity due to hopping between Pr sites, (iii) extrin- 
sic electronic conductivity due to impurities. Compounds 
with the fluorite structure often exhibit ionic conduction 
when vacancies caused by Schottky or Frenkel defects al- 
low the oxygen ions to hop from site to site. However, 
at room temperature the number of thermally induced 
defects is very small, so ionic conductivity in Pr02 is un- 
likely to be important in this temperature range. It has 
been suggested that ionic conductivity in PrO x materi- 
als is significant only above ~ 900Ki2i*22i Therefore, it is 
most likely that the conductivity in Pr02 is electronic. 
Given the small activation energy found here, the mech- 
anism warrants further investigation. 

V. DISCUSSION AND CONCLUSIONS 

We have presented neutron diffraction experiments 
on single crystal Pr02, which reveal an internal distor- 



tion of the oxygen sublattice at Td = 120 ± 2K and 
a related component of the magnetic structure below 
Tn = 13.5 ± 0.2 K. The displacements of the oxygen ions 
are five times as large as those observed in UO^r^ and 
this is reflected in the high value of the transition tem- 
perature. The observation that both the distorted crys- 
tallographic structure and the related component of the 
magnetic structure have a unit cell that is doubled along 
one crystal axis with respect to the fluorite structure sug- 
gests that the displacement of the oxygen ions affects the 
magnetic ordering of the Pr sublattice. 

Our measurement of the specific heat capacity indi- 
cates that the crystal field ground state of the Pr 4/ elec- 
tron is a doublet, consistent with the lowering of the local 
symmetry of the Pr site. This raises the question of what 
mechanism drives the structural distortion. The most 
obvious answer is that the electronic energy is reduced 
by a collective Jahn- Teller distortion at the expense of 
a small penalty in elastic energy. Another possibility is 
that the distortion is a consequence of a quadrupolar or- 
dering of the Pr 4/ orbitals. Intuitively this would seem 
an unlikely mechanism, since the coupling between elec- 
tric quadrupoles would be expected to be too weak to 
cause an orbital ordering at a temperature as high as 
120 K. Nevertheless, there is evidence for a degree of Pr 
4/ — O 2p hybridisation in PrO^r^ 4 ^ an d such a cova- 
lency effect could provide a mechanism for aligning the 
4/ orbitals. 

The splitting of the ground state could also partly ex- 
plain the observation of an ordered magnetic moment 
smaller than that associated with the Ts quartet ground 
state of cubic symmetry. On the other hand, we have 
previously reported a broad continuum of magnetic scat- 
tering in the excitation spectrum of Pr02, which we be- 
lieve arises from a strong magnetoelastic coupling. 8 This 
coupling would tend to quench magnetic degrees of free- 
dom and reduce the size of the ordered moment by the 
dynamic Jahn- Teller mechanism. From the experimen- 
tal evidence available so far it seems likely that both the 
static and dynamic Jahn- Teller effects are important in 
Pr02, but it is not clear which of the two is dominant at 
low temperatures. 

In work that we will report elsewhere^ we have found 
that the application of a magnetic field in the distorted 
phase at T < Td produces striking hysteresis effects in 
both the magnetic structure and the magnetisation. The 
totality of results that have emerged in recent years sug- 
gests that a strong interplay between electronic and lat- 
tice degrees of freedom exists in Pr02 and influences its 
properties, perhaps in unexpected ways. 
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